Discovering Severe Acute Respiratory Syndrome 
Coronavirus 3CL Protease Inhibitors: Virtual 
Screening, Surface Plasmon Resonance, and 
Fluorescence Resonance Energy Transfer Assays 
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An integrated system has been developed for discovering potent inhibitors of severe acute respiratory syndrome coronavirus 
3C-like protease (SARS-CoV 3CL”°) by virtual screening correlating with surface plasmon resonance (SPR) and fluores- 
cence resonance energy transfer (FRET) technologies-based assays. The authors screened 81,287 small molecular com- 
pounds against SPECS database by virtual screening; 256 compounds were subsequently selected for biological evaluation. 
Through SPR technology-based assay, 52 from these 256 compounds were discovered to show binding to SARS-CoV 3CL?". 
The enzymatic inhibition activities of these 52 SARS-CoV 3CL”® binders were further applied to FRET-based assay, and 
IC, values were determined. Based on this integrated assay platform, 8 new SARS-CoV 3CL”” inhibitors were discovered. 
The fact that the obtained IC,, values for the inhibitors are in good accordance with the discovered dissociation equilibrium 
constants (K,s) assayed by SPR implied the reliability of this platform. Our current work is hoped to supply a powerful 
approach in the discovery of potent SARS-CoV 3CL”® inhibitors, and the determined inhibitors could be used as possible 
lead compounds for further research. (Journal of Biomolecular Screening 2006:915-921) 
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T HE LAST OUTBREAK OF THE SEVERE ACUTE RESPIRATORY 
SYNDROME (SARS) epidemic in 2003 has led to thousands 
of lethally infected patients and hundreds of deaths,' and the seri- 
ous aftereffect has caused major medical and economic con- 
cerns. It has been reported that SARS coronavirus (SARS-CoV) 
is responsible for SARS infection,*? and SARS-CoV 3C-like 
proteinase (SARS-CoV 3CL”®) is an attractive target for the dis- 
covery of anti-SARS agents for its functional importance in 
the viral life cycle. To date, varied kinds of SARS-CoV 3CL””° 
inhibitors were discovered,** although no effective SARS-CoV 
3CL”” inhibitor has yet been reported to treat SARS. Therefore, 
it is still a great challenge to explore new chemical classes of 
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SARS-CoV 3CL”® inhibitors that can be possibly used in 
anti-SARS research. 

As has been used in the current drug-discovery process, by 
using the high-performance computation technique to search the 
large chemical compound databases for the identification of pos- 
sible drug candidates, the virtual screening approach is a tech- 
nology that is based on the 3D structure of the target protein. It 
involves the rapid fitting of the chemical library members into 
the active sites of 3D protein structures and is critical to distin- 
guish active from inactive substances at the primary screening 
stage. To date, there have been successful cases of virtual screen- 
ing used in the discovery of lead compound candidates. For 
example, we recently reported that by using virtual screening 
technology with other methods, several new, small molecular 
specific cyclophilin A inhibitors were discovered’; through vir- 
tual screening with surface plasmon resonance (SPR) technology 
against a database containing structural information of more than 
8000 available drugs, Cinanserin (SQ 10,643), a well-characterized 
serotonin antagonist that had undergone preliminary clinical test- 
ing in humans in the 1960s, was discovered to be the SARS-CoV 
3CL”” inhibitor.’ 

Recently, SPR technology has been recognized as a powerful 
tool in monitoring receptor-ligand interactions with advantages of 
no-labeling, real-time, noninvasive measurements and low sample 
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consumption.’ It has been widely used in hit determination,”? 
and more and more successful cases have been published.!°"! 
Fluorescence resonance energy transfer (FRET) technology—based 
assay belongs to another potent system applied in lead compound 
discovery. As a simple and sensitive approach, it has been effec- 
tively used in the detection of proteolytic activity and enzyme 
inhibitor evaluation.*® Usually, FRET-based assay involves an 
intramolecularly quenched fluorogenic substrate that contains 
a donor and an acceptor chromophore. In our previous work, 
we constructed an efficient fluorogenic substrate-based platform 
for screening and evaluation of SARS-CoV 3CL?” inhibitors,’ 
with which several potent SARS-CoV 3CL?° inhibitors were 
determined? and some related enzymatic features of SARS-CoV 
3CL?” were studied." 

In the current work, considering the distinct advantages of 
these 3 above-mentioned methods, and according to the obtained 
successful experiences in enzyme inhibitor discovery,’ we thus 
developed an integrated strategy for exploring novel SARS-CoV 
3CL?” inhibitors. Our work is expected to supply a sound 
approach in the discovery of new potent SARS-CoV 3CL?° 
inhibitors, and the determined inhibitors could be used as possi- 
ble lead compounds for further research. 


MATERIALS AND METHODS 


Materials 


The fluorogenic substrate Dabcyl-KNSTLQSGLRKE-Edans 
was synthesized by Shanghai Sangon Biological Engineering & 
Technology and Service Co. Ltd. (Shanghai, China). The chelat- 
ing column and protein molecular weight marker for sodium 
dodecy] sulfate polyacrylamide gel electrophoresis (SDS-PAGE) 
use were from Amersham Pharmacia Biosciences (Uppsala, 
Sweden). All other chemicals were of reagent grade or ultrapure 
quality and purchased from Sigma (St Louis, MO). Cinanserin 
was synthesized according to the published method.* The tested 
256 compounds from virtual screening were purchased from 
SPECS Company (http://www.specs.net). 


SARS-CoV 3CL”” preparation and activity assay 


SARS-CoV 3CL"° was prepared according to the published 
procedure.'* The enzymatic activity of SARS-CoV 3CL””° was 
measured by cleavage of the fluorogenic substrate Dabcyl- 
KNSTLQSGLRKE-Edans based on the reported protocols.*! 
Briefly, the fluorescence intensity was monitored on a GENios 
microplate reader (Tecan, Mannedorf, Switzerland), and the 
enzyme activity was determined based on the increase in fluo- 
rescence upon continuous monitoring of the reactions in 
96-well black microplates (BMG LABTECH, Offenburg, 
Germany) using wavelengths of 340 nm and 488 nm for exci- 
tation and emission, respectively. 
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Virtual screening 


Virtual screening was performed on a 64-processor SGI 
Origin 3800 supercomputer. The virtual screening protocol 
adopted in this work was similar to that used in the previous pub- 
lication.’ Briefly, the 3D model of SARS-CoV 3CL?” was per- 
formed as the target for screening the SPECS molecular database 
containing 81,287 molecules by a docking approach’* and the 
program DOCK 4.0.'° Residues within a radius of 6 A around the 
catalytic center (His41 and Cys145) were used to construct 
the scoring grids for the docking screen. During the docking 
calculations, Kollman-all-atom charges'’ and Geisterger-Hiickel 
charges'*” were assigned to the protein and the small molecules, 
respectively. Conformational flexibility of the compounds from 
the database was considered in the docking search. In DOCK 
simulation, the ligand-receptor binding energy was approxi- 
mated by the sum of the van der Waals and electrostatic interac- 
tion energies. After initial evaluation of orientation and scoring, 
a grid-based rigid body minimization was carried out for the lig- 
and to locate the nearest local energy minimum within the recep- 
tor binding site. Position and conformation of each docked 
molecule were optimized using the single anchor search and tor- 
sion minimization method of DOCK 4.0. The top 2000 mole- 
cules with the highest DOCK score were further subjected to 
CScore” and Druglike”! calculations. Finally, with visual inspec- 
tion, 256 compounds were selected and purchased. 


SPR technology-based binding assay and kinetic study 


To inspect the possible active compounds with binding affini- 
ties to SARS-3CL?”, SPR technology based Biacore 3000 biosen- 
sor was used (Biacore AB, Uppsala, Sweden). SPR assay was 
performed at 25 °C, and SARS-CoV 3CL”® purified at 3 mg/mL 
was diluted to a final concentration of 25 ug/mL in 10 mM 
sodium acetate buffer (pH 4.3) and immobilized to the CM5 
sensor-chip by using the standard primary amine coupling 
method.” HBS-EP running buffer (10 mM HEPES, 150 mM 
NaCl, 3.4 mM EDTA, and 0.005% [v/v] surfactant P,,, pH 7.4) 
was used as the running buffer during immobilization. The sur- 
face was activated by injecting a 1:1 mixture containing 0.2 M 
1-Ethy]-3-(3-dimethylaminopropyl) carbodiimide hydrochloride 
(EDC) and 50 mM N-hydroxysuccinimide (NHS) for 7 min at 
a flow rate of 20 uL/min. SARS-CoV 3CL”® was injected and 
the surface was then blocked by injecting 1 M ethanolamine at pH 
8.5 for 7 min. The final immobilization level was about 4000 res- 
onance units. 

All the screening assays were performed over the unmodi- 
fied dextran surface and SARS-CoV 3CL”” surface. During the 
SARS-CoV 3CL"° binder assay, the sample concentrations 
were set at | uM and 10 uM. For sample dilution, the pH and 
the concentration of both DMSO and buffer substances in sam- 
ples and running buffer were carefully matched. Each sample 
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assay consisted of a 60-s buffer injection, a 60-s injection, and 
a 120-s undisturbed dissociation phase. The published SARS- 
CoV 3CL”” inhibitor Cinanserin> was used as a positive control 
to confirm the binding activity of SARS-CoV 3CL°, and the 
blank injection was used to check the carryover effects. The 
signal was adjusted for nonspecific binding of the samples to 
dextran matrix by subtracting the signal in the reference chan- 
nel from the signal in the active channel. The binding hits could 
be directly determined by the signal difference between the 2 
set concentrations of the samples. 

To further quantitatively investigate the kinetic binding fea- 
tures of the determined SARS-CoV 3CL”® inhibitors by FRET- 
based assay, SPR technology—based Biacore 3000 was also 
used. During the assay, samples were diluted to 8 different con- 
centrations and injected at a constant flow of 30 uL/min at 
25 °C. Sensorgrams were processed by automatically subtract- 
ing for nonspecific bulk refractive index effects. The kinetic 
parameters were then analyzed using a global data analysis 
program (BIAevaluation 3.1 software, Biacore AB). 


SARS-CoV 3CL?” inhibition assay 


SARS-CoV 3CL?® inhibition assay was performed according 
to the previously described method.’ The inhibition of SARS- 
CoV 3CL™ activity was assayed by preincubating the enzyme 
(1 uM) with the test compounds (50 UM) in 100 mM NaCl, 
20 mM phosphate buffer (pH 7.4) at 4 °C for 2 hours. The reac- 
tion was initiated by the addition of 10 UM of fluorescent sub- 
strate (Dabcyl-KNSTLQSGLRKE-Edans), and the change of the 
fluorescent intensity was monitored on a GENios microplate 
reader (Tecan) using an excitation wavelength of 340 nm and an 
emission wavelength of 488 nm. Control reactions were carried 
out using the same reaction mixture without the test compound 
or the enzyme. Concentration of inhibitor producing 50% inhibi- 
tion (IC,,) calculation was carried out in the absence or presence 
of varying concentrations of the inhibitors by using Origin 7.0 
software (OriginLab Corp., Northampton, MA). 


RESULTS AND DISCUSSION 


Virtual screening and docking results 


The DOCK program was employed as the first step in pre- 
liminarily screening the potential binding molecules against 
SARS-CoV 3CL”” based on the substrate-binding pocket formed 
by residues within a radius of 6 A around the catalytic center 
(His41 and Cys145) of the model structure.”? As a result, 256 
compounds from the SPECS database (containing 81,287 com- 
pounds) were thus selected for subsequent biological assays 
according to their DOCK score, Cscore, Druglike score, and 
finally with visual inspection. The results obviously implied that 
the performed virtual screening has supplied potent compound 
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FIG. 1. Structures of the 8 compounds screened in this study. 


information, although it might suffer from false positives in view 
of some factors like inaccurate scoring functions, crude confor- 
mational statement of target protein and some solvent-related 
terms,’° or false negatives due to the neglect of allosteric effect as 
carried out in our assay. Shown in Figure 1 are the structures 
of the 8 FRET-determined SARS-CoV 3CL”® inhibitors from 
virtual screening. These compounds show good druglikeness 
according to our in-house druglike filter.7! Furthermore, 5 of the 
8 compounds pass the rule of 5 based on the results shown in 
Table 1. 

The docking results (Fig. 2) demonstrated that the discovered 
inhibitors can fit well to the active site of SARS-CoV 3CL?”. 
Figure 3 shows the detailed binding information of DC060251 to 
SARS-CoV 3CL?” as a typical example. Moreover, by using 
DC060251 as a query structure to search MDL Drug Data Report 
(MDDR), Comprehensive Medicinal Chemistry (CMC), and 
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Table 1. Summary of the Properties of the 
8 SARS-CoV 3CL”® Inhibitors 

No. No. No. Rule of 

Molecular Hydrogen Hydrogen 5 (Pass 

Compound Weight Bond Donors Bond Donors CLogP _ or Fail) 
DC060015 558.4 2 a 6.5484 Fail 
DC060087 487.5 1 9 2.9486 Pass 
DC060159 519.4 4 a 5.6938 Fail 
DC060170 541.3 2 9 4.551 Pass 
DC060180 522.1 2 7 4.507 Pass 
DC060245 536.6 1 i 7.056 Fail 
DC060251 470.5 2 7 5.7126 Pass 
DC060256 420.3 1 5 3.2525 Pass 


FIG. 2. The binding model of the 8 inhibitors against SARS-CoV 
3CL” in the active site. DC060251 is rendered in sticks, whereas others 
are shown in lines. 


China Natural Products Database (CNPD), several known drugs 
or preclinical compounds were also identified (Fig. 4), whose 
structures share 90% similarity to that of DC060251. For example, 
compound A (3,3’,3”,3’”-p-Xylene-01.,01,0 ,o’ -tetrayltetrakis(4- 
hydroxy-2H-1-benzopyran-2-one)) from MDDR was a known 
HIV-1 protease inhibitor with IC,, of 1.7 uM and a mean EC,, of 
11.5 uM in an XTT-based anti-HIV assay.”* Compound B (Ethyl 
Biscoumacetate) from CMC was found in 1943 with anticoagulant 
effect, and compound C (Gerberinol 1) from CNPD was the 
main component extracted from a kind of plant gerbera with many 
biological activities.*° These results thus implied that compound 
DC060251 might be of great interest for further study. 


SPR technology-based binder identification 


As a 2nd key process, SARS-CoV 3CL° binder iden- 
tification from the virtual screening results was based on SPR 
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DC060251 
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FIG. 3. Schematic representation of SARS-CoV 3CL”® interaction 
with DC060251. DC060251 and protein residues are rendered in ball- 
and-stick and sticks, respectively. 


technology-based Biacore 3000 biosensor (Biacore AB). 
SPR technology has been widely used to investigate protein- 
macromolecule””** and protein-small molecular compound>’ 
interactions for its obvious advantages such as label-free, sen- 
sitive, real-time, noninvasive measurements, low sample con- 
sumption, and high throughput. Recently, this technology has 
also been applied in the discovery of inhibitors against varied 
enzymes, such as HIV-1 protease,’ human cyclophilin A,’ and 
human 5-lipoxygenase.’° In our previous reports, we ever suc- 
cessfully used SPR technology for discovering SARS-CoV 
3CL?® small molecular inhibitors® and performing SARS-CoV 
3CL?°—telated research." 

In the current work, the SPR-based results showed that of the 
selected 256 compounds from the virtual screening, 52 com- 
pounds were identified as SARS-CoV 3CL”® binders. SPR- 
based technique has supplied a rapid approach in the discovery 
of SARS-CoV 3CL”” binder. However, it is noticed that in some 
cases some hits obtained by SPR assay are probably nonspecific 
binding molecules, which show no inhibitory activities. To solve 
this problem, the screened hits should be further tested by enzy- 
matic inhibitory assays. 


FRET technology-based IC,, determination 


As demonstrated in the SPR-based SARS-CoV 3CL”® binder 
assay, 52 SARS-CoV 3CL?® binders were discovered. To further 
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FIG. 4. Structures of 3 compounds (A, B, C) from MDL Drug Data 
Report (MDDR), Comprehensive Medicinal Chemistry (CMC), and 
China Natural Products Database (CNPD) databases. DC060251 
shows 90% structural similarity to them. 


investigate their possible inhibitory activities against SARS-CoV 
3CL””, the inhibitory effects of the selected 52 binders were 
tested by the FRET method. For this assay, the internally 
quenched fluorogenic substrate Dabcyl-KNSTLQSGLRKE- 
Edans was applied, whose Michaelis constant (Km) against 
SARS-CoV 3CL”” is 49.38 + 5.80 uM (data not shown). This 
value is comparable to those of other recently published fluo- 
rogenic peptides.°”” 

Table 2 lists the yielded IC,,, values of the 8 screened inhibitors 
of SARS-CoV 3CL”°, and Figure 5 shows the dose-response 
curves for 3 typical inhibitors DC060170 (A), DC060180 (B), and 
DC060251 (C). The IC,,, values of these inhibitors are in the range 
of 6.86 to 80.46 UM. 


SPR-based quantitative inhibitor binding affinity 
determination 


It has been known that although the IC. value represents the 
inhibitory activity of an inhibitor to an enzyme, this value is 
subject to the experimental conditions, for example, the struc- 
ture of the substrate and the total concentrations of both enzyme 
and substrate. If any of the experimental conditions change, the 
IC,, value of a given compound will be altered, making IC,, 
from different literatures hard to be compared.*° In comparison 
with IC,,, the equilibrium dissociation constant (K,) for the 
enzyme-inhibitor complex is independent of the experimental 
conditions. SPR-based assay can directly measure K,, values 
between the inhibitors and SARS-CoV 3CL”®. To further inves- 
tigate the K,, values of the inhibitors to SARS-CoV 3CL?®, an 
SPR technology—based Biacore 3000 instrument was also used. 
K,, values of the 8 discovered inhibitors were calculated using 
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Table 2. Summary of Concentration of Inhibitor 
Producing 50% Inhibition (IC) Values and 
Equilibrium Dissociation Constants (K,s) 
for SARS-CoV 3CL?° Inhibitors 


Molecular 

Compound Weight (Da) ICS) (UM) K, (uM) 

DC060015 558.4 49.96 + 4.62 18.60 + 1.18 
DC060087 487.5 80.46 + 2.07 42.50 + 2.10 
DC060159 519.4 42.19 + 3.02 19.14 + 2.03 
DC060170 541.3 6.86 + 0.91 4.23 + 0.81 
DC060180 522.1 9.19 + 2.34 6.50 + 0.23 
DC060245 536.6 21.49 + 2.97 23.91 + 5.63 
DC060251 470.5 9.56 + 1.75 4.78 + 0.29 
DC060256 420.3 56.24 + 3.48 49.91 + 3.98 


the 1:1 Langmuir binding model or the steady-state affinity 
model (BIAevaluation 3.1 software, Biacore AB). Figure 5 
(inset) shows the binding dose-response curves for the 3 com- 
pounds, DC060170 (A), DC060180 (B), and DC060251 (C). As 
indicated in Table 2, the K, values of the compounds range 
from 4.23 to 49.91 uM. As demonstrated in Figure 6, the K,, 
values are well correlated with IC,, values by the correlation 
coefficient (R) of 0.92. 


CONCLUSIONS 


In this work, we developed an integrated approach for dis- 
covering novel SARS-CoV 3CL?® inhibitors. Virtual screening 
method based on the 3D structure of SARS-CoV 3CL"® can 
reliably narrow down the number of potential candidates before 
experimental testing. Although virtual screening might suffer 
from false positives, SPR technology was thus used to deter- 
mine the binders of SARS-CoV 3CL”° and further reduce the 
screening scale and time. Because the hits obtained by SPR- 
based assay may be nonspecific binding compounds, FRET 
method is applied for SARS-CoV 3CL”® inhibitor identifica- 
tion. Our current strategy for discovering SARS-CoV 3CL° 
inhibitor has integrated virtual screening and SPR- and FRET- 
based assays, and with this integrated platform, we have suc- 
cessfully identified 8 novel potent SARS-CoV 3CL”” inhibitors. 
Although the hit enrichment might not always be higher than a 
random selection, the screened active compounds from virtual 
screening could involve more structural information for virtual 
screening based on the compound database that comprises a 
large number of chemical structures, whereas for random 
screening, the obtained active compounds might contain lower 
structural information because of the possible limitation of the 
compound quantity. The obtained results thereby suggested that 
the inhibitor screening is more effective compared with the tradi- 
tional screening when integrated strategy is applied. It is 
believed that such an integrated approach may help speed up the 


www.sbsonline.org 919 


Downloaded from jbx.sagepub.com at TEXAS A & M INTL UNIV on March 16, 2015 


Chen et al. 


Response Una (RU) 
est 88s 


100 150 200 250 300 350 
Time (s) 


Inhibition (%) 
ror) 
° 


10° 10° 10° 10° 10* 10° 
Log(C) M 


100 150 200 250 300 350 
Timne (8) 


Inhibition (%) 
@ 
o 


10° 10° 10° 10° 10* 10° 
Log(C) M 


‘o) 


100 150 200 250 300 350 
Tene (3) 


Inhibition (%) 
o 
o 


10° 10° 10° 10° 10* 10° 
Log(C) M 


FIG. 5. Dose-response inhibition and kinetic assay curves for the 3 
compounds from the screening procedure. IC, values were extracted 
from the curves using nonlinear regression analysis, and dissociation 
equilibrium constant (K,) values were analyzed using a global data 
analysis program of Biacore 3000. A, B, C represent the dose-response 
inhibition curves (Inset: the responding binding curve) of the compounds 
(DC060170, DC060180, DC060251), respectively. 
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FIG. 6. Correlation between the concentration of inhibitor produc- 
ing 50% inhibition (IC,,) and equilibrium dissociation constants (K,) 
for SARS-CoV 3CL”° inhibitors. IC,, values were obtained from 
the fluorescence resonance energy transfer (FRET)—based assay, and 
K, values from biosensor measurements. The symbol (@) represents 
the compounds with correlation. 


drug-discovery process. In addition, it is expected that the deter- 
mined inhibitors could be used as possible lead compounds for 
further research. 
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